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Results of a search for the Flavor-Changing Neutral Current decay B® d — > n + p us- 
ing pp collision data at y/s = 1.96 TeV collected at Fermilab Tevatron collider by the 
CDF and D0 detectors are presented. CDF reports upper limits on B(B^ — » /i + p~) < 
7.5- 1CT 7 and B{B° -> p + p,~) < 1.9- 10~ 7 at the 95% C.L. using 171 pb" 1 . The D0 Col- 
laboration used 240 pb _1 to set an even more stringent limit on the branching ratio for 
B° -> fj,+fj,~ of 5.0 ■ 10~ 7 at the 95% C.L.. 
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1. Introduction 

The Tevatron Run II at Fermilab started to deliver pp collisions at yfs = 1.96 TeV 
in April 2002. The bb production cross section at this center of mass energy is very 
large (rs 150 /ib) when compared to the typical e + e~ cross sections at the T(45)(« 
1 nb) and Z° (w 7 nb) resonances. The bb production rate at the current record 
luminosity at the Tevatron of « 10 32 cm -1 sec -1 is 10 kHz. At hadron colliders there 
is also the advantage of producing all b-flavored species from the light B^ and B® 
mesons to the heavier £?!? and B+ meson as well as b baryons such as A°. On the 
other hand, the problem of the hadronic environment is the high level of background 
due a large inelastic pp cross section of O(50 mb) and the track multiplicity to up 
to 50 tracks/events due to the fragmentation of the hard intersection products plus 
that due to the underlying event or pile up of multiple events. 



2. Theoretical Introduction 

The purely leptonic decay B^ — > is a Flavor-Changing Neutral Current 

(FCNC) process^ I n the Standard Model (SM), this decay is forbidden at the 
tree level and proceeds at a very low rate through higher order diagrams. The SM 
leptonic branching fraction (B) were calculated including QCD correctional The 
latest SM predictionPis, B{B Q S -> = (3-42 ± 0.54) • 10~ 9 , where the error is 

dominated by non-perturbative hadronic uncertainties. The corresponding leptonic 
branching fraction for the B® is suppressed by an additional factor of iVt^/VisI 2 
leading to a SM branching ratio of (1.00 ±0.14) • 10~ 10 . CDF reported in Run I the 
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experimental bound for the branching fraction of (-BS) of B(B° S (B { j) — > < 
2.6 • 10~ 6 (8.6 • 10~ 7 ) at the 95% C.L.B1. 

The decay amplitude of Bd s — > fi + u~ can be significantly enhanced in some ex- 
tensions of the SM. For instance, in the type-II two Higgs Doublet Model (2HDM), 
all contributions from the neutral Higgs sector cancel out and the branching frac- 
tion depends only on the charged Higgs mass Mh+ and tan/3 which defines the 
ratio of the vacuum expectation values of the Higgs field. The amplitude grows 
as tan 4 j3 In the Minimal Supersymmetric Standard Model (MSSM) however, 
B(B® — > tan 6 /3, leading to an enhancement of up to three orders of mag- 

nitudel^Jcompared to the SM, even if MSSM with minimal flavor violation (MFV) 
is considered, i.e., the CKM matrix is the only source of flavor violation. 



3. Event Selection 

For normalization CDF used the b cross section as measured by the experiment in 
Run I; therefore, the preselection criteria were driven by the measurement of the 
cross section with p|? > 6 GeV and \ub\ < 1-0 for events selected by dimuon triggers. 
In addition a set of standard track, muon (p-r > 2 GeV/c, < 0.6) and vertex 
quality cuts were applied. After restricting the invariant mass region of the dimuon 
pair to 4.669 to 5.696 GeV/c 2 , 2940 events survive in a total integrated luminosity 
C of 171 pb" 1 . 

D0 used data collected by dimuon tiggers with two muons of opposite charge 
that form a common secondary 3D-vertex with an invariant mass between 4.5 and 
7.0 GeV/c 2 . Each muon candidate had to have pr > 2.5 GeV/c, |?7| < 2.0 and 
a sufficient number of hits in the central tracking station. To ensure a similar pt 
dependence of the system in the signal and in the normalization channel, ptj?, 

had to be greater than 5 GeV/c. After the preselection, 38,167 events survive in a 
integrated luminosity of 240 pb -1 . 



3.1. Discriminating variables 

Both experiments have chosen a set of similar discriminating variables to best exploit 
the properties of the decay. The long lifetime of the B® meson has been used to 
cut against random combinatoric background due to short lived particles. While 
CDF uses a minimum proper lifetime cr of the B® candidate, D0 is using decay 
length significance Lxy / 5Lxy ■ The fragmentation characteristics of the b quark are 
such that most of its momentum is carried by the b hadron. The number of extra 
tracks near the B candidate therefore tends to be small. The second discriminant 
was therefore an isolation variable, I, of the muon pair, defined as: 

T= Ip(m + aOI m 

track i^B 
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Here, Pi i s the scalar sum over all tracks excluding the muon pair within 

track i^B 

a cone of IS.1Z < 1 around the momentum vector of the muon pair where 

A7Z = \J (A0) 2 + (Ai]) 2 . The CDF definition differs by using the transverse mo- 
mentum of the tracks instead of the momentum. 

As a third variable the angle (A$) between the momentum vector of the muon 
pair and the vector from the primary to the secondary vertex has been used. This 
requirement ensures consistency between the direction of the decay vertex and the 
momentum vector of the B® candidate. 



4. Limit setting procedure and optimization 
CDF 

For a given number of observed events, n, consistent with the background estimate, 
ribg, the upper limit on the branching ratio is determined using: 

B(B s ^^-)< N ^ (2) 

I ■ a Bs ■ a ■ Stotai ■ J Lett 

where N(n, nt, g ) is the number of candidate decays at 90% C.L., estimated using a 
Bayesian approach^! and incorporating the uncertainties into the limit. The inte- 
grated luminosity is £, The B® production cross section is a B o (= j^CTb+^I with 

U =0.100, f u = 0.39^1) and a ■ e total is the total acceptance times efficiency, ob- 
tained from data and Monte Carlo. The factor of two is necessary, since the analysis 
is sensitive to the charge conjugate b hadron. 

The optimization was done with approximately 100 combinations of the dis- 
criminating variables and maximized the a priori expected limit which is given by 
the sum over all possible observations, n, weighted by the corresponding Possion 
probability of the expected n og . The optimization of the discriminating variables 
(cr > 200jUm, A$ < 0.10 rad, 2 > 0.65) predicted 1.05±0.30 background events in 
171 pb- 1 . 

D0 

A random grid search and an optimization procedure^' was used to find the 
optimal cut values of the discriminating variables by maximizing the variable 
P = e^/{a/2 + y/Nb g ). Here, e M 4 is the reconstruction efficiency of the signal 
events relative to the preselection (estimated using MC), and A^g is the expected 
number of background events interpolated from the sidebands. The constant a is 
the number of standard deviations corresponding to the confidence level at which 
the signal hypothesis is tested. This constant a was set to 2.0, corresponding to 
about the 95% C.L. The result (L X y/SL X y > 18.5, A$ < 0.2 rad, 1 > 0.56) of 
the optimization leads to a background prediction of 3.7±1.1 events in 240pb . 

In the absence of an apparent signal a limit on the branching fraction B(B S ) 
can then be computed by normalizing the upper limit of number of events in the 
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B s signal region to the number of reconstructed B^ — > J/tp K events: 

N f 

B(B S ) < -M. • ■ -It^ Bl{B ±) . B 2 (j/^), (3) 
N B± e££ Jb^B u , d 

where iV u i is the upper limit on the number of signal decays estimated from 
the number of observed events and expected background events including sys- 
tematic uncertainties using the Feldman and Cousins ordering scheme for the 
MC integration The number of observed S ± — > J/tpK^ events is Ng± = 
741 ± 31 ± 22, e^/e^ K = 0.247 ± 0.009 ± 0.017 is the efficiency ratio of the 
signal and normalization channels, obtained from Monte Carlo (MC) simulations. 
The fragmentation ratio of a b o r b quark producing a B® and a B^ or B& is 
fb^B a /fb^B n , 4 = 0.270±0.034Elandi3i = B{B ± -> J/tpK^ = (1.00±0.04)-10- 3 
and 2?2 = B(J/tp — > /i/i) = (5.88 ± 0.1)% — are the measured branching ratios of 
the normalization channel. 



5. Results 

Both experiments have used a blind analysis technique, hiding the signal region 
until the analysis was completed. The invariant mass spectra after unblinding are 
shown in Fig.Q] No excess of a signal has been seen in either of the two experiments. 
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Fig. 1. The invariant mass spectrum for the D0 experiment (left side) and the CDF experiment 
(right side) after optimized requirements on the discriminating variables. 



CDF observed one event with an invariant mass of AL+„- = 5.295 GeV/c 2 , thus 
falling into both the B® and B d mass window. Since this is consistent with the 
expected background of 1.1±0.3 events in each of the mass windows (±80 MeV 
around the corresponding PDG mass), 95%(90%) CL limitJ^of B(B° -► fi + fi~) < 
7.5 • 10~ 7 (5.8 • 10~ 7 ) and B{B a d -► < 1-9 ■ 10~ 7 (1.5 • 10~ 7 ) are found. 

D0 found four candidate events in the invariant mass region of ±180 MeV of the 
expectecffSl B° s mass, which is also consistent with the expectation of 3.7±1.1 back- 
ground events. This gives a resulting limil^lon the branching fraction at 95% (90%) 
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CL of B(B° S -> fi + n~) < 5.0 ■ 10" 7 (4.1 ■ 10~ 7 ). This limit assumes that there are 
no contributions from .Bj] — ► decays in the search region. Any non-negligible 

contribution due to B® decays would make the obtained limit on the branching 
fraction smaller. The limit presented for B(B® — > is therefore conservative 

and is even more stringent to constrain models of new physics beyond the SM. 
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